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Abstract. Euler angles determining rotations of a system as a whole are 
conveniently separated in three-particle basis functions. Analytic integration 
of matrix elements over Euler angles is done in a general form. Results for the 
Euler angle integrated matrix elements of a realistic NN interaction are listed. 
The partial wave decomposition of correlated three-body states is considered. 
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1 Introduction 

When expansion methods are applied in three-particle problems matrix elements 
(ME) are to be calculated which are six-dimensional integrals. It is expedient 
to choose Euler angles, representing rotations of a system as a whole, as three 
of the six degrees of freedom and to carry out the integration over Euler an- 
gles analytically. Then one is left with integrations that are three— dimensional 
only. Such low-dimension integrations can very efficiently be done with regular 
quadratures. The elimination of rotational degrees of freedom has been used by 
several groups [jl], [2|, ||] at solving both the three-nucleon Schrodinger dynamic 
equation and that of the method of integral transforms . No restrictions on the 
form of basis functions arise in such an approach, and in conjunction with use 
of an appropriate set of basis functions the techniques are probably the fastest 
and easiest ones to solve most of the nuclear three-body problems of interest. 
Therefore, it seems useful to list the formalism. In Sec. 2 the separation of Euler 
angles is performed for general form three-particle functions. In Sec. 3 analytic 
integration over Euler angles of general type ME is carried out which is applied 
in Sec. 4 to calculate the ME pertaining to a realistic NN interaction. The partial 
wave decomposition of correlated basis functions is considered in Sec. 3 as well. 
Sec. 5 contains comments on symmetry properties of Euler angle integrated ME. 

In Ref. H most of the ME considered in Sec. 4 were also calculated, how- 
ever, our techniques are different. Both the results and techniques of the present 
paper are substantially simpler. Besides, we present the general formula for Eu- 
ler angle integrated ME, and our choice of the body reference frame leads to 
simplifications. 

Contrary to the correlated basis functions case, ME of one-body and two- 
body operators can be calculated in a simple way without the Euler angle inte- 
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gration when e.g. uncorrelated hyperspherical harmonics are used as a basis set. 
However, the present techniques are helpful also in the latter case when ME of 
NNN force or those involving two-body subsystems at solving reaction problems 
are considered. 

While the formulae are written down below for the case of the three-nucleon 
problem they may be applicable also beyond these frames. 

2 Separation of Euler Angles 

We use the Jacobi vectors 

u=-L(r 2 - ri ), v = y| ( r3 - I!±^) (1) 

where r\ are nucleon positions, and we consider three-dimensional rotations 
u,v — > u',v'. We interpret rotations as rotations of a coordinate system so 
that a rotation a — > a' means that {a x ,a y ,a z } = a and {a' x ,a' y ,a' z } = a' are 
coordinates of a vector in the old and new coordinate systems, respectively. Spa- 
tial components Flmi with given total orbital momentum and its projection 
quantum numbers, of three-particle basis functions transform under rotations 
as § 

L 

F LM i(u',v')= d mm>( u ) f lm(u,v) 

M=-L 

where uj denotes collectively three Euler angles parameterizing a rotation and 
Dmm> are ^ ne Wigner D-functions. 

In order to separate the Euler angle dependence of Flm( u , v ) we consider 
rotations u' = u'(u,v), v' = v'(u,v) which are different for different u, v 
vectors, so that u = lo(u,v). We use the inverse relation 

L 

F LM (u,v)= D^ mi (lo)F lm ,(u',v'). (2) 

M'=-L 

At given u and v, we choose the new coordinate system in a way that it cor- 
responds to a body reference frame associated with the u, v plane. The corre- 
sponding body-frame coordinates u'^v^ may be expressed in terms of scalars u, 
v, and (u • v). Treated in this way, Eq. (|2|) becomes an expansion of Flm( u , v ) 
over D-functions depending on Euler angles with coefficients depending on scalar 
variables. The Euler angles determine positions of the body reference frame with 
respect to a laboratory reference frame. 

In particular, a convenient choice j7|] is to direct the z-axis of the body 
reference frame along the vector u, and to place the x-axis of the body reference 
frame in the u, v plane in a way that the projection of the vector v onto this 
axis is positive. Then one has 



u' x = 0, u' y = 0, v! z = u, v' x = vy 1 — t 2 , v' y = 0, v' z = vt (3) 



4 



Elimination of Rotational Degrees of Freedom 



where t = (u ■ v), u and v being the unit vectors in the directions of u and v. 
As a result, one comes to the expansion 

L 

F LM (u,v)= Yl DMM>(v) FLM> (u,v,t) (4) 

M'=-L 

where the Euler angles oj parametrize the above defined rotation into the body 
reference frame, and the notation Flm_(u, v, t) stands for the quantities 
Flm{u\v') = FLAt(u' x ,u' y ,u' z ,v' x ,v' y ,v' z ) taken at the space point from @: 

Flm (u, v, t) = F LM (u' x = 0,u' y = 0, u' z = u ,v' x = vy/l - t 2 , v' y = 0,v' z = vt). 

(5) 

Writing down the formulae we shall consider F^m to be real which is commented 
below in connection with Eq. (|i~9[). 

The coordinates u> = {a, (3, 7}, u, v, t ranging between their natural bound- 
aries realize a mapping onto the whole phase space (u, v). The correspondence 
between the points in the phase space and the a, (3, 7, u, v, t coordinates is a one- 
to-one correspondence everywhere except for the exceptional sub-areas where 
the u, v plane is undefined and hence the body reference frame cannot be de- 
fined. These are sub-areas where vectors u and v are collinear, or u = 0, or 
v = 0. Many sets of Euler angles correspond to the same point in the phase 
space if the point belongs to these sub-areas. For the case of integrands which 
are non-singular in the exceptional sub-areas the contribution of these sub-areas 
to ME vanishes (not speaking of the fact that the phase space volume element 
suppresses the contributions of u = and v = sub-areas, see Eq. (0) below). 
Therefore, one can transform the integration over u, v into the integration over 
a, P, 7, u, v, t and apply Eq. 

Our basis functions Flm are regular single-valued functions of the points 
of the phase space. This is in contrast to D-functions themselves which do not 
possess such property in the above mentioned exceptional sub-areas. While this 
does not cause problems when Eq. (||) is used for integration purposes, use of 
products of D-functions and functions of internal variables as basis functions 
would be a disad vantage. 

Let us also comment on the aspect of symmetry with respect to particle per- 
mutations. A well-known way to get totally antisymmetric three-nucleon basis 

1 In addition, it is easy to see that Eq. (^) remains valid in the exceptional sub-areas where 
u and v are collinear, and where v = 0. Both sides of Eq. (Q) turn to Ylm{/3, a)f(u, v) in 
these sub-areas, and integration over them using the coordinates a, /3, 7, u, v, t leads to correct 
results. 

2 At given L and M values D-functions form a complete set of 2L + 1 functions of three coordi- 
nates determining the position of a system as a whole. Constructing the basis set Flm one, for 
example, can use another set of 2L + 1 functions instead of D-functions which have the same 
property but which are basically polynomials of the components of the vectors u,v and are 
regular in the whole phase space. These functions are to be multiplied by functions of three 
internal variables. This basis proved to be very efficient in atomic physics 
When such a basis is used to calculate ME of some simple operators there is no need to apply 
Eq. (^). The integrands in the left-hand side of Eq. (|^) below can be decomposed in a closed 
form into the sum of components with given L values. Only the L — component contributes 
to the right-hand side of (W), and no summation over M" is then required, see tel. 
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functions is to combine spatial basis functions that are components of irreducible 
representations of the three-particle permutation group with spin-isospin func- 
tions of conjugated symmetry. The corresponding spatial basis functions are often 
constructed via application of symmetrization operators to complete sets of ini- 
tial "simple" functions. It is convenient, and sufficient, to use the expansion (||) 
of symmetrized functions involving one and the same body reference frame de- 
scribed above. In contrast to the basis functions themselves, this reference frame 
does not possess symmetry properties with respect to permutations of three par- 
ticles. The arising symmetrized functions F aym are expressed in terms of initial 
functions F with the help of relations of the type 

F sym (u, v)=J2 C P F(Pu, Pv) (6) 
p 

where P are particle permutations, and Cp are given coefficients. Hence, the 
required body-system components (|5|) of symmetrized functions are obtained 
with the help of such type relations applied at the points (vf, v') given by Eq. 
(^). One thus needs to form linear combinations of the corresponding vectors 

Pu' = apu' + bpv', Pv' = cpu' + dpv' 

and to calculate initial functions at the points (Pu' , Pv'). When, alternatively, 
the basis functions antisymmetric with respect to nucleon permutations are ob- 
tained via a direct antisymmetrization of angular momentum coupled products 
of space functions and spin-isospin functions the procedure is similar. 

ME between functions of the form (|j) will be calculated below. The volume 
element is 

dr = dudv = dujdTi n t, dr^t = u 2 v 2 dudvdt. (7) 

Here the relationship dudv = dojdt, see Appendix, is used. The overlap-type 
integrals are 

J dTG* L , M ,d U)V)t F L M = S LL i5 M m'8tt 2 (2L + 

r L 

x / dT int G L m" (u, t)O u ,v,t F L M" (u, v, t). (8) 

J M"=-L 

Here O u<Vt t is an operator acting on u,v,t, and the orthogonality properties 
( A.37Q are used. 

3 Matrix Elements of Tensor Operators 

To calculate ME it is expedient to decompose the operators into sums of tensors 
irreducible with respect to rotations. Performing integration over Euler angles 
analytically is then possible in practically important cases. The following rela- 
tionships are useful to this aim. 

1. Let Ok K be the irreducible tensor operator || of the rank k. We obtain its 
Euler angle integrated reduced ME (Gx/||Ofc||.Fi) defined as usual, 

{Gv M '\O kK \F LM ) = CfX(2i' + iy l/2 (G L/ \\6 k \\F L ), 
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where CYm^ k are the Clebsh-Gordan coefficients. The result is 

(Gv\\6 k \\F L ) =^ 2 (2L' + iy l l 2 J dr int Y, ClmL Gl'M' 6 kK F LM . (9) 



MM'k 
M'=M+k 



Here O kn F^M means the quantity O^kFlm taken at the space point from (||). 
The relation (J9|) is useful when one is able to calculate OukFlm in a simple 
way. The cases when the operator O kn includes multiplications by functions and 
differentiations fall into this category, cf. the next section. 
To get Eq. (||) we define the states 

(O*0&fc = £ C^ K 6 kK F LM . (10) 

k+M=M 

Since O kK is an irreducible tensor these states possess definite rotational quan- 
tum numbers Lq, Mq. One therefore can expand them, as in Eq. (|j), in terms 
of products of D-functions and the body-system functions {OF)^ Mo {u, v, t) de- 
fined similar to Eq. (||). Writing down the relation inverse to ( ]T0| ) and utilizing 
this expansion gives 

o kK F LM = c^d^; m ,(u) (of)^ M) . 

L M M> 

Use of this relation, the D-function expansion of Gl'M', and Eq. (||) gives 

dTG* L , M ,d kK F LM = 87r 2 (2L' + l) _1 Cfj$j. K dT int £ G L i M " {OF) k jj M n . 
J J M"=-L> 

Substitution of the body-system version of Eq. flu]) here, 

{Q F ) k VM" = Y C LMk K OkKFLM_, (11) 
M+k=M" 

gives Eq. (§). 

More complicated ME can also be integrated over Euler angles. Let us com- 
ment on e.g. the ME of the square of a three-nucleon Hamiltonian H 2 = (T+V) 2 
entering the well-known Temple lower bound for binding energy. ME of T 2 are 
calculated like the overlap integrals (||) with Flm —> TFlm , Glm —> TGlm , 
ME of TV + VT are calculated using the same replacements and Eq. (||), and 
ME of V 2 are given by the sum of contributions of the form 

((G L2 xs 2 )j M I {0 k2 ■ Z k2 ) (6 kl ■ Z kl ) | (F Ll <p Sl )j M ) ( 12 ) 

where O k and E k are the rank k = 0, 1, and 2 tensor operators, {O k ■ £ k ) is their 
scalar product ||, and (Fl (g) <Ps)jm symbolizes the LSJ coupling. The opera- 
tors O k and U k act in the coordinate subspace and the spin-isospin subspace, 
respectively, and they are irreducible with respect to rotations. 

To calculate ME (||) one performs recouplings, see Appendix, leading to 
spatial and spin states of the type {OF)ffl , (OG)%\ {^) S S % and (27x)|£ 
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where the notation is similar to that of Eq. (10). Subsequent application of Eqs. 
d) and (|) to the states {OF)^\ (OG) 1 ^ 2 gives the following expression for 
the ME (ID, 



Y I(M 1 ,k u M 2 ,k 2 )J(M 1 ,k u M 2 ,k 2 ), (13) 

I(Mi,Ki,M 2 ,K 2 ) = J dT int d k2K _ 2 GL 2 M 2 Qfc^^LiMi , 

8vr 2(_ 1 )5 1 -5 2 ^ (2(S+1) | S 
LS { 1 



S L J 11 S L J 

Si h | | L 2 S 2 k 2 



J(M 1 ,K 1 ,M 2 ,K 2 ) 

1 

Here M s is arbitrary. 

2. The partial wave decomposition 

F LM (U,V) =J2yLM(u,v)f hh L(u,v), (14) 
ilia 

Y^(u,v) = C^ lhm2 Y hmi (u)Y hm2 (v) (15) 

mi+m2=M 

of correlated three-body states could also help to calculate ME between them. 
This applies e.g. to ME of a non-local NN interaction, or NN interaction with 
/-dependent form factors. To perform the decomposition we need to calculate the 
overlap integrals 

f khL (u,v) = J dudvYlif(u,v)F LM (u,v). (16) 



If used directly, Eq. (|l_y) implies the four-dimensional integration to get / at 
each u, v point. But separating Euler angles and performing integration over 
them analytically one can replace the four-dimensional integration with a single- 
dimensional one. Indeed, similar to Eq. (|8|), using the Eq. (Q)-type expansion for 
both factors in the right-hand side of Eq. (16) we obtain that 

f hhL {u,v) = ^ 2 {2L + l)- 1 f 1 dt Y Y^(t)F LM (u,v,t), (17) 

J 1 M=—X 

where Y^(t) are the body-system values of Y^*(u,v), 

Yjjm = A l t l ^Y hM {cose = t,<p = 0), A h L % = [(2h + 1)/(4^)] 1 / 2 C^ M , 

(18) 

and A = min(/2, L). 

If in the expansion (14) the functions /; 1 / 2 l(^, v) are realf^] then the body- 



system functions Flm are real as well. This is seen from the body-system version 



of Eq. (14) 



F LM (u,v,t) = Y Y lM(t)fhhUu,v) (19) 
hh 



3 In the usual case of T-invariant Hamiltonians the expansion coefficients are real for such a 
choice of basis functions. 
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since, with the choice made of the body reference frame, the functions 



entering (|19| ) are also real. 

4 Matrix Elements of NN Interaction and Kinetic Energy 

The basis functions are antisymmetric with respect to particle permutations so 
it is sufficient to calculate the ME of an interaction between a pair of particles. 
We consider here such an interaction, e.g. V(12), that includes the operators 

I 2 , k 2 , (l-(s 1 + s 2 )), S 12 = 3(n-a 1 )(n-a 2 )-((T 1 -cT 2 ), 
Sf 2 = 3(1 • o-i)(i • tr 2 ) - P(<7i • * 2 ), = 3(fe • • <r 2 ) - k\a l ■ er 2 ) (20) 

times factors V[r) and isospin and spin projection operators. Here n = r/r, 
r = r 2 — ri, k = (l/2){k 2 — k\) = —id/dr, I = r x k, and <Tj are the Pauli 
matrices. Most frequently used versions of the realistic NN interaction are of this 
form. We shall use Eq. (||) to obtain ME of the above operators integrated over 
Euler angles. Initial ME can alternatively be calculated using Eqs. ([lj]), ( |l~7l) but 
the first of these procedures is faster. 

Within e.g. the approximation of equal proton and neutron mass (denoted 
with m) ME of kinetic energy are expressed in terms of those of the k operator. 
Indeed, kinetic energy T of A nucleons in the center-of-mass system can be 
written as (2mA)~ l J2i<j(fii ~Pj) 2 i where pi are the nucleon momenta. Then 
ME of T can obviously be calculated as 

(G\f\F) = (A- l)(4m)- 1 (G|(ft - p 2 f\F) = (A - l)(h 2 /m)(G\k 2 \F). 

Final expressions for the Euler angle integrated ME of the operators (^) are 
given by Eqs. @, fl23|), @, J27|), @, @, and @ below. ME of the Is force 



and of further operators in (|20| ) are expressed in the usual way H in terms of 
the reduced ME of spatial components of those operators. In the Is force case, 
the corresponding Euler angle integrated reduced ME {GLi\\V{r)i\\Fj}j of the 

operator V{r)l^ are required. Equation @ gives 



Gu\\V{r)l\\F L )=& 



3tt 2 (2L' + 1)~V2 f dr mt V(r) ^ C^ lt G UM > IjJlm - 

MM', n=±l 



M=±l 
M'=M+im 

(21) 



Here the quantities I^Flm are defined as above. One has IqFlm = which is 
taken into account in (|2lD, For /j, = ±1 these quantities are of the form 




'-^--T2{aK ±i ^^ FLM (22) 

where the derivatives are taken at the space point from (||). The derivatives can 
be calculated numerically, and high accuracy is attained in this way. To perform 
this for symmetrized basis functions one needs to apply Eq. @ at space points 
(it, v) that are obtained from Eq. (||) by the replacement of u' x = 0, or u' y = 0, 
with small non-zero values of these quantities. 
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Only real parts are to be collected in the right-hand side of (22) when the 



functions f\ x i 2 L are real in Eq. (|l~4|). Indeed, the quantities 1±\Flm are real in 

/■hfa 
LM ■ 



this case which is seen using the expansion of 1±\Flm over y' 1 ' 2 



When calculating ME of the I 2 operator it is convenient to express them in 
terms of the quantities that already appear in ME (|2l|), (The same is done below 
also in Eqs. ([24] ) and (^).) For this purpose one writes down the integrand in 
the form V(r) Y.^{^G L m)*{1 I iF L m), expresses Z M F LM and l^G LM in terms of the 
Eq. ( |lO|) type quantities, uses subsequently Eqs. (||) and (||), and uses expressions 
of Eq. (11) type and the orthogonality property of Clebsch-Gordan coefficients. 



This gives 

J dTG* L , M ,V(r)l 2 F LM = 

$LL' $ MAI' 

8^ 2 (2L + 1)- 1 / dT mt V{r) Y hG LM „ iyF LM „ . (23) 

J M",fi=±l 

The Euler angle integrated ME of the k 2 operator are obtained as 

J dTG* L , M ,k 2 F LM = (1/2) J dTG* L , M ,A u F LM 

Sll'Smw f , ( dG* LM dF LM P \ Slu^mm'^ 2 

= 2 J dT l^T^T + G ^ Flm ) = 2L + 1 

x f dT int {Y f [(d/du) G LM » }\(d/du) F L M» }+u- 2 Y l » G LM" l»F LM » }.(2A) 

J M" M",n=±l 

When ME of the operator l/2[V(r)k 2 + k 2 V(r)\ are considered the integrand in 
the last line of is obviously to be replaced with 

V(r){Y i [(d/du) G LM » ][(d/du) F LM , / ] + u- 2 Y l ^°LM" l»F LM „ ] 
M" M",n=±l 

+2- 1 / 2 V'(r)Y{ G LM ,, [(d/du)F LM „] + [(d/du) G LM » ] F LM „ }. 

M" 

The quantities Svz, S 1 ^, and 5jf in ( p0| ) can be written as 

Y x l?4? = (3/2) Y 

i,j=l p=-2 
(2) (2) 

where X^- and E^- are the Cartesian components of a space tensor and spin ten- 

(2) (2) 

sor, respectively, and X^ and E^ are the corresponding spherical components. 
The latter are linear combinations [0] of the Cartesian ones, and in particular 
X { 2) = xi 2 J, E^ ] = Ef). In our case 

(2) 

E^' = (l/2)(<Tii<T 2 j +aij<j 2 i) - (l/3)S ij (er 1 • <r 2 ). 
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One thus has 



= I E A^5(-ir, 

s l A = | E ^E^i-ir, s% = § E Kp^-iy. 

H=-2 n=-2 



The tensors , £^ 2 \ and /C^ are 

Mf = A(n/5)^ 2 Y 2 ,(n), 
= i(7r/5) 1 / 2 {y 2 ,(i)} sym , 4 2 ) = 4(vr/5) 1 / 2 ^(fc), (25) 

where y 2ll (x) = x 2 Y 2 u(x) are the solid spherical harmonics. The \i = compo- 
nents of these tensors are 

A^ 2) = 3n 2 - 1, £ {2) =3P Z -P, 4 2) =3fc 2 -P. (26) 

To perform integration over Euler angles in the reduced ME of the first of 
the tensors (p5|) one uses Eq. (||) . One has A/"j? = 25^o that gives 



(^||^(r)Ar( 2 )||Gr) = 16vr 2 (2L / + l)- 1 / 2 / rfr, nt y(r)^C % Glm - 

(27) 

The Euler angle integrated reduced ME of the two other tensors (|25|) will be 
obtained via calculating the corresponding ME of the operators (p6|). Proceed- 
ing in the same way as at the derivation of Eq. ( p3| ) we come to the following 
expressions 



(G L2 \\V(r)C^\\F Ll \ 

8vr 2 / V(r) dr int £ f(L 1: L 2 , M x , M 2 , m, fi 2 ) l^Gi* M 2 IjmFLiMi , (28) 

J w =±l,Afi,M2=±l,Af 2 
fii+Mi=H2+M2 

(G L2 ||(l/2)[y(r)/C( 2 ) +K^V{r)]\\F Ll ) = 
8^ 2 / dr itrt ]T f{L 1 ,L 2 ,M 1 ,M 2 

,Ml>M2) |^(^) V^2^2M 2 V^i-^LiMi 

jUiAfi^a-Wa 
Ml +Mi=/i 2 +M 2 

+2- 1 /V(r) (^ i0 V M2 Gl 2 a/ 2 FlxMx +^q G L2 m 2 V^^A/j ] (29) 
which contain the quantity 

(2 L 2 + 1)V 2 

-<L 2 M 
y LiM20 



/(Li,L 2 ,M 1 ,M 2 ,/x 1 , At2 ) = ^ L2M 5(Li,L2,Mi,M 2 ,M,/ii,/x 2 ), (30) 



llX^fnr' i i\-l/^Z/Af r»L'Af .oZ/M' riL'M' \ 
6 y2-^y ZL " ~r l ) ( - / LiAflO°L 2 A/10 L 'LiMil A ti ( - y L 2 M 2 l/i 2 ^ 



S{Li,L 2 , Mi,M 2 ,M, HiiH 2 ) 

' • (31) 



2Li+l 
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It is implied here that \L\ — L 2 | < 2. The value of M in ( pOj ) is arbitrary provided 
that C^ 20 / 0. In © M' = /zi + Mi = /i 2 + M 2 . In © V = d/dr. The 
quantities V u,Flm and V^Glm are real at the same conditions as above. 



The summation in (|31|) can be done in a closed form. For this purpose, let 
us perform a transformation of the sum. It follows from the general structure 
of the ME § that the equality S = C^ m S' with S' being independent of 

M holds true. Then multiplying S by C L ^M20> summing over M and using the 
relation J2m (c_LiM2o) = (^2 + -0/5 one concludes that the quantity (|30| ) may 
be represented as 5(2L 2 + 1) ^ 2 <Sb, where So = J2m ^LiM20^- ^- n d due ^° ^ ne 
relation X)m ^lIm20 = the ^ as ^ term from ( |3l| ) is eliminated in the course of 
the summation. As a result, one obtains 

f(Li, L 2 , Mi, M 2 , /xi, ^ 2 ) = 15(2L 2 + 1)-^ £ ^^O^X^X,*. ( 32 ) 

L'M' 

g(L 1 ,L 2 ,L') = (21/ + 1) 1 C^^ 20 C^ 1 % w C^ 2 m 10 .(33) 

M 

(The summation over M' in (^) is a formal one.) Now we perform the summation 
in (|33|) with the help of Eq. 8.7.(13) in f| obtaining 



We substitute this into ( p2[ ) and perform the summations with the help of 
Eq. 8.7.(32) in f§. Finally, we obtain that in ©, (p£ 



f(Li,L 2 , Mi, M 2 , /ii, /x 2 ) = (e^^C-ir+^C^^i^Cjg;^. (34) 

It is implied here that /ii + Mi = /x 2 + M 2 , otherwise / = 0. 

The formulae in this section have been verified in three-nucleon computa- 
tions. 



5 Symmetry relations 

We write down the formulae for (components of) states possessing definite parity 
(— 1) ■ The relation 

Fl,—m = {-1) I+L+M Flm (35) 
holds true. According to ( |35| ) FL t M=o vanishes when I + L is odd. To obtain ( |35| ) 
one notes that in Eq. (|l9j) Y^_ 2 M = (-l) h+h+L+M Y^ . 

Let us assume that the quantities O^kFlm have definite parity (—1)* (that 
is related to the parity of Flm and that of the operator) . Then 

Ok,-KF L ,- M = (-iy +L+M+k+K Ok K F LM . (36) 

Eq. ( |36[ ) is obtained expressing its left-hand side in terms of the functions 
(OF)^ M from Eq. ([ll]) and applying Eq. (|35|) to them. 
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The integrands in Eqs. (||), fll7|), fl27D, and ( p3[ ) may be written, up to the 
notation of summation, as 

J2 X ( M ) and Y, X ( M > ti- 
nt Mu 

Due to Eqs. ( |35| ) and (|^) one has 

X(-M) = (-1) I+I 'X(M), X(-M,-v) = (-1) I+I 'X(M,V), 

where (— 1)* and (— 1)^ are the parities of the two states entering an ME. Then, 
when parities are the same, one can apply e.g. the relations 

^I(M) = I(0) + 2^X(M), £X(M, M ) = £X(0, M )+2 £ X{M,n) 

M M>0 Mfi A* M>0,fi 

to speed up the summations. When parities are opposite the sums vanish as it 
should be. 

The integrands in Eqs. (|l3|), (|2l|), (pi), and (p9|) may be written as 



X] X(M,M',fi) and £ X(M U fn, M 2 , » 2 ). 

MM'n=±l M lfJ , 1 M 2 H2 
M'=M+fl M\+n= M 2+l J -2 

One has, with the same definitions of (— l) 1 and (— l) 1 ' as above, 

X(-M,-M',- M ) = X(M, M', fi), 

X(—Mi,—fix, —M 2 , — £i 2 ) = (-1) 7+/ 'X(M 1 ,/X 1 ,M 2 , M2 ). 

When one of the summation variables above is expressed in terms of the others 
this leads to the sums of the form 

53 Y(m 1 ,m 2 ), Y(m 1 ,m 2 ,m 3 ) 

with the similar property Y(—mi,...,—mi) = (— l) I+I 'Y(mi, . . . , raj). Then, 
when parities are the same, at a particular choice of summation variables one 
can apply e.g. the relations of the type 

Y / Y(M, f i)=J2y(0,^) + 2 ]T Y(M,fj,), 

Mfj, /u M>0,(U 

^ y(M 1 , m , / u 2 )= ^ y(0, m ,/i 2 ) + 2 53 F(M 1)Ml ,// 2 ). 

Ml/ii/12 Ml ^2 Mi>0,/ii//2 
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Appendix 

The known relationship dudv = dujdt is obtained as follows. Let us recall the definition, see 
, of the Euler angles uj = {a,/3,~/}: the angles a and /3 are the azimuthal angle and polar 
angle determining the direction of the vector u, and the angle 7 corresponds to a rotation of an 
intermediate reference frame around its z-axis coinciding with the u direction. (This rotation 
brings the z-axis of the intermediate reference frame into the u,v plane.) Let us perform the 
integration over dv choosing some fixed reference frame with the z-axis directed along u. Then 
one can write dudv = da sin /3d/3dtdip where tp is the azimuthal angle of the vector v in the 
chosen reference frame. But, by the definition of 7, ip differs from 7 only by a constant, so that 
dip = dj. This gives the desired expression taking into account that j(| duo = da sin ftdfid'y. 
The D-function orthogonality relations we use are || 

/8tt 2 
dujD MzM> 2 ^) D M 1 M[ H = + l &LxL^M 1 M 2 & M [M' 2 ■ ( A -37) 

To obtain we need to expand the states of the form 

[Ok ■ Ek) I (Fi ® ips) JM ) 

over those of the form 

\((OF) l *®(Zv)t) JM ). 

The expansion coefficients are the same as if Ok and Sk were not operators but states of 
independent subsystems with angular momenta k. Hence, these coefficients are the recoupling 
coefficients for four angular momenta ||, i.e. basically 9j'-symbols. Our initial states include 
scalar products so that two of the four angular momenta are coupled to zero angular momentum. 
Therefore, the 9j-symbols reduce to 6j'-symbols in our case. 
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